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Abstract. We apply the optimized effective potential method (OPM) to the 
multiplet energies of the 3d™ transition metal atoms, where the orbital dependence of 
the energy functional with respect to orbital wave function is the single-configuration 
HF form. We find that the calculated OPM exchange potential can be represented by 
the following two forms. Firstly, the difference between OPM exchange potentials of 
the multiplet states can be approximated by the linear combination of the potentials 
derived from the Slater integrals F' 2 (3d, 3d) and F 4 (3d, 3d) for the average energy of 
the configuration. Secondly, the OPM exchange potential can be expressed as the 
linear combination of the OPM exchange potentials of the single determinants. 
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1. Introduction 

The density functional theory (DFT) has a important problem with the calculation 
of transition metal complexes. In spite of the recent development of the exchange- 
correlation functionals, it is difficult to evaluate the multiplet energies using the Kohn- 
Sham method, because these electronic structures are not always written in a single 
Slater determinant. Several methods have been proposed to calculate the multiplet 
energies using DFT. One of these methods has been proposed by Ziegler 126 , Wood [21] 
and von Barth |25j . They have carried out the computation of the multiplet energies 
using the diagonal sum rule, in which the sum of the multiplet energies is equal to 
the corresponding sum of the single determinant energies. The method reproduces the 
multiplet energies of p electron system, but the method does not give correct multiplet 
energies for the d electron system. [1] 

The time dependent density functional theory (TDDFT) has been recently 
developed to calculate excited states including the multiplet states. The method has 
been applied to p electron systems [22] within the adiabatic approximation. The 
non-adiabatic TDDFT calculation, in which the exchange-correlation kernel becomes 
frequency dependent, is still difficult. 

Another method was formally discussed in the Gorling's paper [2E] in which the 
energy functional depends on the multiplet states. They applied the optimized effective 
potential method to the calculation of the multiplet energies of p electron atoms. 14 
The optimized effective potential method (OPM) was firstly proposed by Sharp and 
HortonpQ. The method was applied to the Neon and Carbon atoms by Talman et. 
al.[2]. Gorling pointed out that OPM is formally regarded as DFT because the orbital 
dependent energy functional implicitly depends on the density [21] • The OPM potential 
can be evaluated by the singular integral equation called as the OPM equation. Talman 
et. al.[21E] and Engel et. al.|22] solved the one-dimensional OPM equation for an atom, 
in which the angular components of the OPM equation for an atom are integrated out 
using the spherical symmetry. To efficiently calculate the OPM potentials for atoms 
and molecules, Gorling et. al. expanded the OPM potential using some base functions. 
[IH E3 CHI [21] The base expansion method gives the same total energy as calculated by 
Talman's method. However, in the base expansion method, the OPM potential shows 
spurious oscillations. To solve the problem, some methods have been proposed but still 
are not applied to d electron atoms. [TBI E3 120] 

An early work of OPM for the multiplet states is the OPM calculation by Aashamar 
[H] who applied OPM to the ground multiplet of the atom of which atomic number is less 
than 85. The OPM potentials for the p electron ground and excited state multiplets are 
calculated by Nagy [BJ using the KLI method. [HJ E] Furthermore, Aashamar also applied 
OPM to the multiplet energies of p electron atoms using the multi-configuration HF 
form. [3] The total energies of the multiplet states have been discussed in these studies, 
but the potentials of the multiplet states are not presented. Sala and Gorling have 
been developed the LHF method to approximately calculate the exchange only OPM 
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potential. [TTJ H21 H3| Recently, the LHF method is extended to the open shell atoms 
and molecules. ^U] They applied the LHF method to the ground and excited multiplets 
of p electron atoms and molecules. However, in these papers, they did not treat the 
excited multiplet energies of the d electron atoms, which are difficult to calculate using 
the DFT. Furthermore, these studies give incorrect result caused by the approximation 
methods such as the KLI or the LHF method. 

In the present paper, we apply OPM to the ground and excited state multiplets 
of the 3d n electron atoms using the total energy functional of the single-configuration 
HF [31 form. Following Talman and Engel[2*| I22j. we solve the one-dimensional OPM 
equation correctly. The solution is considered to give the most reliable result among 
all the available methods. We show that the differences between the OPM potentials 
of the multiplet states can be approximated by the linear combination of the potentials 
derived from F 2 (3d, 3d),F 4 (3d, 3d) Slater integrals. We numerically show that the OPM 
exchange potential of multiplet state is approximated by the linear combination of the 
OPM exchange potentials of the single determinants. 



2. Optimized Effective Potential Method 

In the present paper, we calculate the OPM potential for eigenstates of an transition 
metal atom. Each state is characterized by a definite value of the total orbital angular 
momentum L and of the total spin angular momentum S. These multiplet energies 
can be expressed as E LS,S , where s is the seniority number. The multiplet energy 
level is degenerate for possible directions of total momentums L and S. Then, each 
multiplet wavefunction ^(LSL Z S Z , s) can be written as a linear combination of single 
Slater determinant functions <&{L z S z ,v) which are not always eigenf unctions of the 
atom. L z and S z are z component of the total momentums L and S, respectively, and 
v specifies different states with the same L z and S z component. Each total energy of 
the single determinant is represented as E * * ,v . 

Furthermore, we consider "average energy of the configuration" E av defined by 
Slater. [321 EHj For the 3d n multiplets, _E av is expressed as 

BT = (1 ° ~ Q ^ )!n! ^2(2L + 1)(2S+ l)E LS ' s . (1) 

LS,s 

To generally express these three kind of energies, 

E LS,s^ E L z S z ,u and ^ 

we 

introduce where £ stands for {LS, s},{L z S z , u} and {av}. The total energy 
is divided into the 1-body part E\ and the 2-body part i?f; 

E* = E* + El (2) 

The 1-body part is constructed from the kinetic energy of electrons and the electron- 
nucleus Coulomb energy. 

Et = y %2q i I*(n i li), (3) 
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J e (n i Z i ) = J P^r^P^dr, (4) 

i_l 2dr 2 + 2r 2 r J ' 1 } 

where represents the number of the electrons in the orbital P„.j.(r). Z is the atomic 
number. rii and U are the principal quantum number and the angular quantum number, 
respectively. In this study, we suppose that the OPM potential V r Q PM (r) is spherically 
symmetric. The single electron orbital is determined by the following equation 

(£i + Vi PM (r))Pi k (r) = 4pt ik (r), (6) 

where ef is the single electron energy. P„./.( r ) * s the product of the distance from the 
nucleus r and the radial component of the single electron wave function. 

The 2-body part E\ is composed of the electron-electron interaction energy 
expressed as 

E l = \ E 4^') + \ E b UGl(ij), (7) 

ijk i7^j,k 

where the coefficients a\- k and b\- k depend on the multiplet states, the single 
determinants or the average energy of the configuration. The Slater integrals F^{ij) 
and G^(ij) are defined as 

F*(ij) = JdrJ dr'pt ik (r)Pl k (r) ^P^P^r'), (8) 

G\{ij) = JdrJ dr'Pijr)P^(r)^P^y)P^ h (r'), (9) 

where r< is the smaller of r , r' and r> is the larger. 

The OPM potential VQ PM (r) is determined by requiring that be minimized for 
all Pn ik obtained from equation (|Bj). This results in 

\ / " AV ' = o (10) 

where 



£ ps ^ jy^ "»»* v - SP^ t (r 



Using equations and the variation of the normalization condition 

5P € , (r) 
<^ PM (r') 



! P r H u( r ? dr = 1 with respect to Vg PM (r'), that is, f P^ k {r) ^^. dr = 0, we can 



rearrange equation (fTUj) as 



V I dr' ( 2 qi V^ M {r)Pi l (/) - ~44t1 = °- < 12 > 
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The functional derivative — , appeared in equation (ITUl) is calculated using the 

^opmW 

perturbation theory: 

dV OPMv) jjki e j e i 

The functional derivative f^ 2 in equation (fT2|) is calculated as 
5El 



J 2 



5Pl 



Uik V J jk jk 

i 



E 2 Ww^4(0 +E4aW)<,>') ) ( 14 ) 



where X^(ij,r) is defined as X%(ij,r) = J Pn i i i ( r ')Pn i u( r ')^^^ r ' '■ Substituting 
equations (fT3j) and (fUjl into equation (fT2j) . we obtain 



= EE / ^(v.O^(OGf(»',O^S i i 4 (r)cir'. (15) 

Equation (fT3j) is called as the OPM equation. The OPM potential can be obtained by 
self-consistently solving the OPM equation and the single electron equation (JHJ). 

3. Method of Calculations 

To solve the OPM equation (fT5|). we modified the code developed by Fischer et. al. [31J. 
For the calculation of the Green's function, we use the following expression 

GHr,r>) = Pi^Qi^) - P^(r')<,(r) - P n ^(r)<,(r') 

+ Cl ik Pi ili {r)Pi ik {r'l (16) 



t dP K (r) 

where Q nil .{r) is the second solution of equation (0) satisfying — ^ — Q n j 



^Pi h (r) = 1, and 



Jo 

/oo 
(17) 

with C iih = I PiiS r ^US r ) dr - To eval uate P^ h (r) and <2^.(r), we solve the single 
electron equation (JBJ using the Numerov's method with the transformation h = log(Zr) 
and P ni ii(h) = Pn^ir) / y/r . The lower bound of h is set to -6.0 and the upper 
bound, which depends on the kind of atom, is set around 7. The step size Ah is 
0.005. The integration in the OPM equation is approximated as the sum of the 
product of integrand and step size. Using this simple approximation with the step 
size Ah = 0.01, the OPM equation is approximated to the set of linear equations for 
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the OPM potential. The OPM equation (|T5)l determines the solution VQ PM (r) only up 
to a constant. The constant is fixed by the physical requirement lim^oo Vq PM (t) = 0. 
In our program, we fix the constant using the asymptotic form of the OPM potential 

- Efe I( a NNk + b NNk)^T P n N l N ( r ') rfr ' at r = 9 aU ' Where N iS tlle illdeX fOT H0M0 - 

The self-consistent procedure converges satisfactorily in about 12 iterations if the average 
of the initial and final VopmM is taken at each iteration. The SCF iteration has 
converged when max^y/g^AP^J < 1.0 x 10~ 8 y/ZN occ , where iV occ is the number of 
the occupied orbitals. [3T] 

In OPM, the single electron energy of HOMO Cqpm * s ec L ua l to the HF single 
electron expectation value for HOMO e^p jHJ. In our results of the calculations for the 3d 
transition metal atoms, the differences between the single particle energies (I^opm - 6 hfI) 
are less than 0.005 au. 



4. Results and Discussions 

The OPM exchange potential Vj^(r) is defined by subtracting spherical average of the 
Coulomb potential from the OPM potential Vq 



Vi(r) = V ? PM (r) - / ^dr>, 



OPM 



r] 



r> 

where the radial density ffi(r) is defined as 

pHr) = }Z^ P U r ) 2 - (19) 

i 

Firstly, we apply OPM to the average energy P av defined by equation (0). The 
average OPM potential Vop M (r) and average single electron orbital P^.( r ) are obtained 
by the self-consistent solution of the OPM equation (|15[) and the single electron equation 
(Pi. The average radial density p av (r) is calculated using equation (fT9|). In DFT, 
the exchange energy is expressed as E DFTex \n(r)] where n{r) is an electron density. 
To calculate averaged DFT potential, we used spherical averaged density p av (r)/47rr 2 
denoted as n av (r). Then, the average energy for DFT is defined as 

^dft = J2 + / ^f^p- drdr' + E DFTex [n av (r)]. (20) 

i J \ \ 

The average density is calculated by the self-consistent solution of the Kohn-Sham 
equation with spherical exchange potential SEdf s ^}^ — ■ 

In the bottom part of figure [H we show the Xa exchange potential (a = 2/3), the 
Becke's GGA (B88) exchange potential and the OPM exchange potential V^(r) of 
Mn 2+ . The OPM exchange potential is close to the Xa and B88 exchange potential. The 
B88 exchange potential diverges at the nucleus, whereas the OPM exchange potential 
does not diverge. 

The radial density p av (r) is shown in the top part of figure ^ There are some 
kinks in the OPM exchange potential. The positions of the kinks correspond with the 
positions of troughs in the radial density. The kinks in the B88 exchange potential are 
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smoother than that in the OPM exchange potential. For the Xa potential, the kinks 
are not appeared. 

In table [TJ we show the average total energies (Equation (JTJ)) using HF method 
and OPM. The average total energy of DFT (Equation (jlOJ)) using Xa and Becke's 
GGA (B88) exchange functional is also tabulated in table HJ The total energy of the 
OPM method is more closer to the HF total energy than that of the B88 and Xa. 
The difference between the HF and OPM energies is order of 0.005 au throughout the 
transition metals. 

Secondly, we discuss the OPM exchange potentials for the 3d n multiplet states. 
The 2-body part of the total energy of an atom (equation can be represented as the 
linear combination of F^(ij) and G^(ij). For the 3d n multiplet states, the coefficients 
a 3d,3d,2 an d a 3d,3d,4 depend on the multiplet states, while the other coefficients do not. 
Therefore, the 2-body part of the total energy is expressed as 

El = El + 4 dj3d)2 Ff (3d, 3d) + 4 i3dj4 F c 4 (3d, 3d), (21) 

where Eq is the 3d n multiplet independent part of 2-body energy. 

Since E%, Fj?(ij) and Eq are functionals of density n£(r) = p^(r)/Airr 2 , we can 

define the potential for JSf, Fj?(ij) and Eq as V| 2 (r) = ^jfc, V^(r) = ^iff and 

Vj (r) = j^qyj, respectively. In the following, we derive the equation to determine 
these potentials. We define £^[-P„./.] as a general expression of El, Fg(ij) and Eq. The 
corresponding potential is defined as = s ^f, r \ ■ From equation (JHJ), the single electron 
orbital P^. is considered to be a functional of Vqpm- Furthermore, we regard Vq PM as 
a functional of n£. Following chain rule of functional derivative, we get 



e ' J J ^ SPi , <<■') 6VL„,(r") * 



QPMV r// ) ^22) 



<5n'(r) 

.\jn i izjZju uy i txi 

' v J J 8\ 

hand side over r and r" . we obtain 



Multiplying the both sides of equation fl2H) by J dr—^-^—, and integrating the right- 



/ drY 2 qi P* (r) 6P ^ kir) Vkr) = [ dr'V ^^^^{23) 
where we used the following relation: f dr—p-^— = f dr V. 2fljPf, (r) - r ih \ j- . If we 

° V OPMV ) "OPM^ > 

substitute £| for ^[P^.] in equation Hfity. we get the OPM equation (JT2J). Furthermore, 
if we substitute F^(ij) and Eq for S^P^,], equation (J23J) defines the potential V^ k and 

Vj o , respectively. Since equation (j2HJ) defines Vg(r) only up to a constant, we set the 
boundary condition linv-^x, Vq PM (t) = 0. Equation (J23|l is solved using the method 
described in sectional 

Functional derivative of equation ()21j) with respect to n LS ' s (r) leads the 
decomposition of the OPM potential for the multiplet state V^pj^(r); 
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In the top part of figure 121 we show the radial density p L5 > s (r) of each multiplet states 
of V 2+ . In the bottom part of figure [21 we display the difference of the radial density 
p LS ' s (r) — p av (r). Comparing the top and the bottom parts of figure we find that 
the difference between the radial densities of the multiplet states are much smaller 
than the radial densities. So we neglect the difference of the radial components P^, v 
between the average configuration and the multiplet states, and we substitute (r) 

^3d,3d 

and V% (r) for Vjg' 9 (r) and VT S (r) in equation (|24)l . respectively. The difference 

3d, 3d ^3d,3d ^3d,3d 

between the OPM exchange potential for the multiplet states V e ^ s ' s (r) and that for the 
average energy V^(r) can be approximated as a linear combination of V^i (r) and 
V%\ (r): 



3d, 3d 



F 4 

3d, 3d 



V^> s (r) - VZ(r) (a 3 L d 5 3 s dj2 - a? )V% Jr 

' ' 7 5 3d, 3d 

LS,s _ av - 

3d, 3d 



+ (<3 S d,4-«W^ M «- (25) 

' ' 3d, 3d 

Neglecting the multiplet dependence of the orbital, we can replace the orbital P^, t . in 
equation (jlj) and (|2~Tj) to P™i- Then, the 1-body and 2-body energies are approximated 

as 

In order to discuss the equality of equations (|23|) and (|2l))l. we compare the 
approximate OPM exchange potentials of V 2+ multiplets with exact ones. The dotted 
line of figure El exhibits the difference between the OPM exchange potential of the 
multiplet states and that of the average energy calculated from equation (|25jl . In the 
solid line of figure El we show the same difference calculated from exact equation . 
The dotted line of figure El close to the solid line. Small differences come from the neglect 
of the multiplet dependence of the radial wave functions. 

In third column of table El we show the approximate total energy of V 2+ multiplets 
calculated from equation ()2f)j1 . For comparison, the exact OPM total energies are 
shown in second column of table 121 If the total multiplet energies are deviate from 
the average of configuration, the approximate equation (f2T)j) poorly predicts the total 
energies. However, the differences between the approximate and the exact total energies 
are no more than 0.002 au. For the other 3d n transition metal atoms, the differences 
between the approximate and the exact energies are no more than 0.005 au. Therefore, 
we can conclude that the equation (|25|) and (|26|) work as the approximation formula of 
the exchange potential and the total energy for the multiplet states, respectively. 

In stead of showing the OPM potentials for many multiplet states, we illustrated 
Vp2 (r) and V^J (r) multiplied by the occupation number g 3 d in the top parts 

-^3d,3d ^3d,3d 

of figure and figure respectively. As atomic number increases, the potentials 
Q3dVp2 (r) and qz&VSi (r) shrink and become high. This corresponds with the spread 

-^d^d ^3d,3d 

tendency of the electron density or the single electron orbital. In the bottom parts 
of figure |U and figure El we show functional derivative of Ff v (3d, 3d) and F^ v (3d, 3d) 
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with respect to P^J(r), respectively. Maximum of ^ i?a gp^' 3d ' > is approximate to that of 
Q3dVpl (r), but the shapes of the potentials are different. The two peaks are observed in 

^3(1, 3d 

Q3dVp2 (r), whereas, a single peak in 5FsjV S pt' 3d ^ ■ Furthermore, at the nucleus, 5Fav S pt' 3d ^ 

3d, 3d 3d 3d 

and q3dV%2 (r) are zero and nonzero value, respectively. 

-^3d,3d 

Finally, we discuss the relation between the OPM exchange potential for the 
multiplet state and that for the single determinant. The multiplet states ^(LSL Z S Z , s) 
are represented as the linear combination of the single determinants Q(L Z S Z , v). If the 
radial components of the wavefunctions do not depend on the multiplet states, the single 
determinant §{L z S z ,v) is represented as 

${L M S M ,u) = Y,Y. a LSL z s^{LSL z S Z)S ). (27) 

LS s 

The 2-body energy of the single determinant Eo zSz)V is the expectation value of V\ ■ — 
with respect to $>(L Z S Z , v>), which can be expressed as 

E LzS z ,u = ^^2 a V LSL z S z ,s a LSL z S z ,s< E 2 S,S i ( 28 ) 
LS s,s> 

where £'^" s ' s = f ^/*(LSL Z S Z , s')J2ij j~^{LSL z S z , s)dri ■ ■ ■ oItn, and r is the spatial 
and spin coordinates. Selecting some combination of single determinants, we solve 
simultaneous equations (jZEj) with respect to the 2-body part of the multiplet state. The 
solution of the simultaneous equations is approximately valid because the orbital wave 
functions depend on multiplet states or single determinants. 

E FS ' S ~ Pl z s z ,»E^ s -\ (29) 

L z S z ,u 

Replaceing £\ . y- appearing in the derivation of (f2*9"|) with total Hamiltonian, we obtain 
approximate equation for the total energy. 

e ls,s „ f3 LzSz ,uE L ^ v . (30) 

L z S z ,u 

From the variation of equation (J2*9")l with respect to n LS ' s (r), the OPM potential 
for the multiplet state V^p^r) = ^fjL^ is expressed as the linear combination of the 

OPM potentials for the single determinants V p M z ' (r) = j^psjns '■ 

V$Z(r) *> Y ^s^Vifu^r). (31) 

L z S z ,f 

Following equation ([31jh we evaluate the approximate OPM exchange potential for the 
multiplet state. Figure El presents the difference between the OPM potential for V 2+ 
multiplet state and that for the average energy. The dotted line exhibits the approximate 
difference potential calculated from equation (|31|). We select the single determinant 
appearing in equation (f3Tj) as |2, 1, — 2|, |2,1,1|, |2,1,0|, |2, 1,— 1|, |2, — 2,2|, |1,0,1|, 
|2, 1,— 2|, 1 2 , 1 , 1 , |2,1,2|, where the number represents the z component of orbital 
angular momentum. The absence and presence of the line over the number stand 
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for spin up and down state, respectively. In the solid line of figure |H1 we show the 
difference potential calculated from exact equation ([24)1 . We find, from figure El that 
the approximate difference potential ( dotted line ) is very close to the exact one ( solid 
line ). 

In the fourth column of table El we show the approximate total energy calculated 
from equation (jSOjl . We compare the approximate total energy with exact OPM total 
energy presented in second column of table 121 The difference between the approximate 
total energy and the exact one is no more than 0.004 au. For the other 3d transition 
metal atoms, the difference between two energies is no more than 0.01 au. Therefore, 
the exchange potential for the multiplet state can be evaluated using the theory such 
as the density functional theory, where the electronic structure is based on the single 
determinant, 

5. Conclusion 

We applied OPM to the multiplet energies of the 3d transition metal atoms using the 
total energy functional of the single-configuration HF |3T] form. The calculated OPM 
exchange potential can be approximated by the Xa and B88 exchange potentials in 
r > 0.1 au. For near nucleus region (r < 0.1 au ), the OPM exchange potential 
strongly deviates from the Xa and B88 exchange potentials. The difference between 
the OPM exchange potential for the multiplet states and that for the average energy 
is quite smaller than the OPM exchange potential. The OPM exchange potential for 
the multiplet state of the 3d transition metal atom can be represented as the linear 
combination of Vi 2 , V^ 4 and VL. We find that the OPM exchange potential 

3d, 3d ^3d,3d 

can be approximated by the linear combination of multiplet independent potentials 
Vp2 and VpJ . This result indicates that the total energy functional can be 

^3d,3d ^3d,3d 

represented as linear combination of the terms which are the product of the multiplet 
independent quantities derived from the Slater integrals (Ff v (3d, 3d) and F^ v (3d, 3d)) 
and the multiplet dependent quantities (%^ d2 , CL^ d ^ dA ). 

As the representative of the potentials for the many multiplet states of the 3d 
transition metal atoms, we discussed the features of the potentials VS? and VSl 

3d, 3d ^3d,3d 

We find that the potentials VS£ and V^i shrink, and become high as the atomic 

^3d,3d *3d,3d 

number increases, which is the same trend as observed in the wave function of the 3d 
transition metals. 

The multiplet state of the 3d transition metal atom is represented as the linear 
combination of the single determinants. From our calculation result, we find that 
the OPM exchange potential of the multiplet state can be approximated as the linear 
combination of the OPM exchange potentials of the single determinants. The result 
might be similar to the methods of Ziegler [2E] , Wood [21] and von Barth |2S] • However 
our result shows that the d electron multiplet state energy should be evaluated by 
including multiplet dependency to the total energy functional. 
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Table 1. The average energy of the configuration (in au ) for 3d™ multiplets calculated 
by Hartree-Fock (£g v F ), OPM (^ V PM ), Becke 88 {E^ 88 ), and Xa (E^ a ,a = 2/3) 
method. 







rpscv Jpscv 
-^OEP — -^HF 


rpscv Vpscv 
^B88 ~ -^HF 


jpscv Jpscv 
— -^HF 


Ti 2 + 


-847.6927 


0.0061 


0.0096 


2.9047 


V 2+ 


-942.0952 


0.0063 


-0.0059 


3.0478 


Cr 2+ 


-1042.4387 


0.0065 


-0.0257 


3.1876 


Mn 2+ 


-1148.8609 


0.0066 


-0.0492 


3.3252 


Fe 2+ 


-1261.4995 


0.0067 


-0.0754 


3.4613 


Co 2 + 


-1380.4916 


0.0068 


-0.1034 


3.5968 


Ni 2 + 


-1505.9743 


0.0068 


-0.1320 


3.7330 



Table 2. The total energies (in au ) of V 2+ multiplet states calculated by the HF 
method, OPM, equation (|2(j[) and equation (|30|l . 



Multiplet 


HF 


OPM 


Eq. (JUJ) 


Eq. (jUD 


4 F 


-942.1799 


-942.1733 


-942.1727 


-942.1733 


4p 


-942.1128 


-942.1065 


-942.1065 


-942.1064 


2 H 


-942.0902 


-942.0840 


-942.0840 


-942.0840 


2 G 


-942.1124 


-942.1061 


-942.1060 


-942.1060 


2p 


-942.0245 


-942.0181 


-942.0177 


-942.0161 


2p 


-942.0902 


-942.0840 


-942.0840 


-942.0842 


?D 


-941.9693 


-941.9628 


-941.9612 


-941.9606 


3 D 


-942.0507 


-942.0444 


-942.0442 


-942.0478 
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Figure 1. The radial density p av (r) of Mn 2+ calculated using the Xa, B88 and OPM 
exchange potentials for the average energy (top). The Xa, B88 and OPM exchange 
potentials for the average energy of Mn 2+ (bottom) . 
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Figure 2. The radial densities p LS ' s (r) of V 2+ multiplet states (top). The difference 
between the radial density of the multiplet state and that of the average energy: 
pLS,s( r ) - p ^{r) (bottom). 
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Figure 3. The difference between the OPM exchange potential of V 2+ multiplet state 
(Kx S ' s ) and tnat of tl:i e average energy (V™): V^ 5 ' 6 ' - V™ '. The dotted line is the 
difference potential calculated from equation Q25[l. The solid line is the same difference 
calculated from exact equation (|24H . 
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Figure 4. The potential of the Slater integral F^ v (3d, 3d) multiplied by the occupation 
number: VSS (top). The functional derivative of the Slater integral i 7 ? (3d, 3d) 

-''ad, 3d 

with respect to SF »^ 3d) (bottom). 
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Figure 5. The potential of the Slater integral F^ v (3d, 3d) multiplied by the occupation 
number: q^dVp^ (top). The functional derivative of the Slater integral F^ v (3d, 3d) 

-^3d,3d 

with respect to P^: <5f ^g a d v ' 3d) (bottom). 
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Figure 6. The difference between the OPM exchange potential of the V 2+ multiplet 
state (Vj^ s - S ) and that of the average energy (V£): Kx S,s - KT- Tne dotted line 
is the difference potential calculated from equation l|31l) . The solid line is the same 
difference calculated from exact equation 



